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We explore the potential of measuring the solar neutrino oscillation parameters in the 
proposed gadolinium loaded Super-Kamiokande (SK-Gd) detector. Gadolinium dissolved in 
_C ■ water can detect neutrons much more efficiently than pure water. This imparts the detector 

the ability to observe electron type antineutrinos, transforming Super-Kamiokande into a 
huge reactor antineutrino detector with an event rate approximately 43 times higher than 
that observed in KamLAND. We simulate the reactor antineutrino data expected in this 
high statistics detector. We use these prospective data to study the precision with which 
the solar neutrino oscillation parameters, Attiq and sin 2 Q , can be determined i) with the 
SK-Gd detector, and ii) by combining the SK-Gd data with the global data on solar neutrino 
oscillations. For comparison and completeness the allowed regions of Atoq and sin 2 00, 
expected to be obtained from the data of the solar neutrino and KamLAND experiments, 
are also presented. We find that the SK-Gd experiment could provide one of the most precise 
(if not the most precise) determination of the solar neutrino oscillation parameters Am and 
sin 2 00 . 
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1 Introduction 



At present we have compelling evidences for oscillations of both solar (z/ e ) and atmospheric (y^ 
and Pfj) neutrinos, driven by nonzero neutrino masses and neutrino mixing [1, 2]. They have 
been accumulated over a long period of time in the experiments with solar neutrinos [3, 4, 5, 
6, 7, 8, 9, 10], Homestake, Kamiokande, SAGE, GALLEX/GNO, Super-Kamiokande and SNO, 
in the KamLAND experiment with reactor antineutrinos [11], and in the studies of the fluxes 
of atmospheric neutrinos by the Super-Kamiokande collaboration [12, 13, 14]. Indications for 
oscillations of neutrinos have also been obtained in the K2K long base-line experiment [15]. 

Evidences for oscillations of the solar v e , induced by nonzero neutrino masses and neutrino 
mixing [2], have been reported first by the pioneering Davis et al. (Homestake) experiment [16, 17]. 
They have been confirmed and reinforced later by Kamiokande, SAGE, GALLEX/ GNO and Super- 
Kamiokande experiments. The evidences for mixing and oscillations of the solar v e have been made 
compelling during the last three years by the data from the SNO solar neutrino and KamLAND 
reactor antineutrino experiments [8, 9, 10, 11] (see also, e.g., [18]). Under the assumption of CPT- 
invariance, the observed disappearance of reactor v e in the KamLAND experiment, in particular, 
confirmed the interpretation of the solar neutrino data in terms of u e — > u^ T oscillations, induced by 
nonzero neutrino masses and nontrivial neutrino mixing. The KamLAND results [11] practically 
established the large mixing angle (LMA) MSW solution (see, e.g., [8]) as unique solution of 
the solar neutrino problem. The solar neutrino and KamLAND data, including the salt phase 
SNO results, strongly favor the low-LMA solution with [10, 19, 20] Am^ ~7x 10~ 5 eV 2 and 
sin 2 9 12 ~ 0.3, Am^ = Atoq and 6> 12 = 9 Q being the neutrino mass squared difference and mixing 
angle driving the solar neutrino oscillations. The high-LMA solution (see, e.g., [11]), characterized 
by Am| x ~ 1.5 x 10~ 4 eV 2 and similar value of sin 2 12 , is only allowed at 99.14% C.L. by the data 
[19]- 

Strong evidences for oscillations of the atmospheric (z/^) have been obtained in the Super- 
Kamiokande experiment from the observed i) Zenith angle dependence of the multi-GeV and 
sub-GeV /x-like events [12, 13], and ii) the recently reported observation of an "oscillation dip" in 
the L/E— dependence of the (essentially multi-GeV) \i— like atmospheric neutrino events [14] 1 , L 
and E being the distance traveled by neutrinos and the neutrino energy. As is well known, the SK 
atmospheric neutrino data is best described in terms of dominant two-neutrino — > v T [p^ — > u T ) 
vacuum oscillations with [13, 14] Am 2 tm = Am| x ~ (2.0 — 3.0) x 10 -3 eV 2 and maximal mixing, 
sin 2 2# atm = sin 2 2^ 2 3 ~ 1.0. The observed dip is predicted due to the oscillatory dependence on 
L/E of the — > v T and — > v T oscillation probabilities. This beautiful result represents the 
first ever observation of a direct effect of the oscillatory dependence on L/E of the probability of 
neutrino oscillations in vacuum. 

Information on the third mixing angle #13, present in the case of 3- neutrino mixing which 
is required to describe the solar and atmospheric neutrino (and KamLAND) data in terms of 
neutrino oscillations, is provided essentially by the short baseline reactor antineutrino experiments 
CHOOZ and Palo Verde [21]. Using the range of allowed values of Am 2 tm found in the latest SK 

1 The sample used in the analysis of the L/E dependence consists of /i— like events for which the relative 
uncertainty in the experimental determination of the L/E ratio does not exceed 70%. 
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data analysis [13, 22], a combined 3-u oscillation analysis of the solar neutrino, CHOOZ and 
KamLAND data gives (at 99.73% C.L.) sin 2 13 < 0.074 [19]. A somewhat more stringent limit 
was obtained in the global analysis of all the available solar, atmospheric and reactor neutrino 
data, performed in [23]: sin 2 13 < 0.050 (99.73% C.L.). 

Neutrino flavor oscillations have also been claimed to have been observed in the LSND exper- 
iment [24]. This interpretation of the LSND results is being currently tested in the MiniBOONE 
experiment [25]. 

After the remarkable progress made in the last few years in establishing the existence of neu- 
trino oscillations, one of the main goals of the future experimental studies of neutrino mixing is to 
measure with high precision the parameters which drive the solar and atmospheric neutrino oscil- 
lations, Am 2 ,!, Am 31 , sin 2 0i2 and sin 2 023- The potential of the current and future experiments for 
high precision determination of the solar neutrino oscillation parameters, Am^ and sin 2 0i 2 , was 
studied recently in [26, 27, 28, 29] (see also [30]). The precision in the measurement of Am^ and 
sin 2 0i2, which can be reached in the Borexino and the so-called LowNu experiments, was explored 
in [26, 27]. The impact of the prospective increase in statistics of the KamLAND data on the 
determination of Am^ and sin 2 12 was investigated in detail in [19, 26, 29]. If the e + — spectrum 
measured at KamLAND is simulated at a point in the low-LMA region, the allowed 3a area in the 
high-LMA zone reduces in size in the case of 0.41 kTy of data, and disappears if the statistics is 
increased to 1.0 kTy [19]. In this case Am^ will be determined with high precision. If, however, 
the spectrum observed in KamLAND conforms to a point in the high-LMA region, the conflicting 
trend of solar and KamLAND data would make the high-LMA solution reappear at 90% C.L. 
and the determination of Ara^ would remain ambiguous [19]. In both cases sin 2 12 cannot be 
determined with high precision. Actually, reaching high accuracy in the measurement of sin 2 i2 
(e.g., 10% error at la) is a rather challenging problem. It can be solved by performing a dedicated 
reactor v e experiment [26, 28], or by the future solar neutrino experiments LENS and XMASS, 
aiming to measure the pp neutrino flux [31, 32]. 

Recently it was proposed in [33] to dope the SK detector with Gadolinium by dissolving 0.2% 
gadolinium trichloride in the water 2 . The added gadolinium would make it possible to detect the 
neutrons, released in the v e capture on protons, with a relatively high efficiency This would allow 
the SK experiment, in particular, to detect the reactor antineutrinos coming from the numerous 
powerful nuclear reactors located in Japan, thus transforming SK into a huge reactor antineutrino 
detector with an event rate approximately 43 times higher than that observed in KamLAND. 
In this paper we study the prospects of measuring the solar neutrino oscillation parameters by 
observing the reactor antineutrino oscillations in the proposed gadolinium loaded SK (SK-Gd) 
experiment. Since Am^ <C Ara^ and the mixing angle 0i 3 is restricted to be relatively small, 
the third (heaviest) neutrino with definite mass is expected to have negligible impact on the 
determination of the values of the solar neutrino oscillation parameters 3 in the experiment of 
interest and we perform our study in the framework of the 2z/— mixing scenario (sin 2 13 = 0). 

We begin in Section 2 with an overview of the currently allowed ranges of values of the solar 

2 The authors of [33] called "GADZOOKS!" the resulting detector. We will use in what follows the abbreviation 
SK-Gd for it. 

3 The effect of sin 2 6*13 can become important if sin 2 20 12 is measured with a ~ 10% precision, which, as we will 
see, is higher than the precision which can be reached in the SK-Gd experiment. 
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Figure 1: The 90%, 95%, 99%, 99.73% C.L. allowed regions in the Am 2 21 - sin 2 12 plane. The 
left-hand and right-hand panels show the allowed regions obtained in x 2 — analysis respectively of 
the global data from the solar neutrino experiments and of the solar neutrino and KamLAND 
data. 



neutrino oscillation parameters by the existing global solar neutrino and KamLAND data. We 
present further the regions in the Am^ — sin 2 (9 12 plane, which are expected to be allowed after tak- 
ing into account future higher statistics and lower systematic data from SNO and from KamLAND 
experiments. The SK-Gd detector is considered in Section 3, where we discuss the assumptions 
made to simulate the data in this proposed modification of the SK experiment, as well as the 
procedure used to statistically analyze the simulated data. We investigate further the potential of 
the SK-Gd experiment in reducing the uncertainties in the values of Am 21 and sin 2 # 12 . In Section 
4 the expected results from a combined analysis of the solar neutrino data and the prospective 
SK-Gd simulated spectrum data are given. In Sections 3 and 4 we also discuss the bounds on 
Am 21 and sin 2 # 12 one would obtain if some of the reactors in Japan would be "switched off". 
Section 5 contains the main conclusions of our study 



2 Solar Neutrino Oscillation Parameters from the Solar 
Neutrino and KamLAND Data 

We begin by reviewing the current status of determination of the solar neutrino oscillation pa- 
rameters. We present in Fig. 1 the regions in the Am 21 — sin 2 # 12 plane, allowed by the global 
solar neutrino data (left-hand panel), and by the combined data from the KamLAND and solar 
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Table 1: The range of parameter values allowed at 99% CL. and their corresponding spread. 



5 



1 .7e-04 



1.3e-04 



9.0e-05 



5.0e-05 



1 .Oe-05 



Global Solar with lower 
SNO CC and NC errors 



Global Solar with lower 
SNO CC and NC errors 
+ 1.3 kTy KamLAND 
data simulated at 
low-LMA 



Global Solar with lower 
SNO CC and NC errors 
+ 1.3 kTy KamLAND data 
— simulated at high— LMA 





0.2 



0.3 0.4 

sin 0,, 



0.2 



0.3 0.4 
sin"0., 



0.2 



0.3 



0.4 



0.5 



sin 2 9. 



Figure 2: The 90%, 95%, 99%, 99.73% C.L. allowed regions in the Am^ - sin 2 12 plane. The left- 
panel shows the the areas obtained in a x 2 analysis of the global solar neutrino data with assumed 
errors of 5% and 6% on the values of the CC and NC event rates respectively, measured in the 
SNO experiment. The right panel shows the allowed region expected if 1.3 kTy KamLAND data 
were added to the prospective global solar neutrino data (see text for details). 



neutrino experiments (right-hand panel). The global solar neutrino data used in the analysis 
include the total event rates measured in the Homestake [3] and SAGE+GALLEX+GNO [5, 6] 
(combined) experiments, the Super-Kamiokande 44 bin Zenith angle spectrum data [7], the 34 bin 
day-night spectrum data from the D 2 phase of the SNO experiment [9] and the charged current 
(CC), neutral current (NC) and elastic scattering (ES) data from the salt phase of SNO [10] 4 . 
In what concerns the KamLAND results, we use the 13 binned spectrum data released by the 
KamLAND collaboration [11] 5 . Fi gure 1 was obtained using the updated standard solar model 
(BP2004) results on the solar neutrino fluxes and the associated errors [36]. The range of allowed 
values of the parameters Am^ and sin 2 #12 are shown in Table 1. Also shown are the % spread, 
which is defined as, 

1 dmax Q"min /-, \ 

spread = , (1) 

®"max 0"min 

where a max (a m i n ) are the maximal (minimal) allowed value of the parameter a at 99% C.L. 

The SNO detector can simultaneously measure the CC and NC rates induced by the flux of 
solar S B neutrinos having energy E v ft 6.5 MeV and E v ft 2.2 MeV, respectively. While CC rate 

4 We refer the reader to [18, 19, 34] for further details of the solar neutrino data analysis. 

5 For details of the statistical analysis procedure used for handling the KamLAND spectral data see ref. [35]. 
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depends both on the 8 B flux normalization and the solar neutrino survival probability, the NC rate 
is determined, in the case of flavor oscillation involving only active neutrinos, solely by the 8 B flux 
normalization. This gives SNO the ability to determine the average v e survival probability and 
hence the oscillation parameters with a relatively good precision. In particular, the combination 
of CC and NC rates in SNO can be effectively used to obtain stringent upper limits on the allowed 
values of both Am^ and sin 2 9 i2 [37, 38]. The upper limits on Am| x and sin 2 9 i2 follow from the 
fact that the ratio of the CC and NC event rates observed in SNO is significantly smaller than 
0.50 [10]. The smaller the ratio of the CC and NC event rates observed in SNO, the smaller the 
maximal allowed values of both Am^ and sin 2 # 12 [37, 38]. For values of Am^ and sin 2 # 12 in the 
low-LMA region, the matter-enhanced transitions of 8 B neutrinos inside the Sun are adiabatic and 
the S B v e survival probability, which affects the CC event rate observed in SNO during day-time, 
is given approximately by P ee sin 2 #12. Thus, improvement in the precision with which the CC 
and NC rates are measured in SNO would lead to a diminishing of the maximal allowed values 
of Amjj and sin 2 #12, provided the mean values of these two observables do not change. This 
is illustrated in the left-hand panel of Fig. 2 where we show the regions of values of Am^ and 
sin 2 #12, which would be allowed by the global solar neutrino data if the experimental errors in the 
CC and NC event rates measured by SNO are reduced respectively to 5% and 6% 6 , while the 
mean values of the two rates coincide with those found in the salt phase of the SNO experiment 
[10]. We note that the maximal allowed values of both Am 2 x and sin 2 # 12 would be smaller, while 
the minimal values essentially do not change. Thus, the uncertainties in the values of Am^ and 
sin 2 612 would be reduced with respect to the currently existing ones. 

Much better precision in the determination of the value of Am^ might be achieved in the 
KamLAND experiment. In Table 1 we present the range of allowed values of this parameter 
and the corresponding spread expected if 1 kTy data of KamLAND is combined with the current 
solar neutrino data (see refs. [30, 40] for details of the KamLAND prospective data analyses). The 
expected precision in the determination of Am^ from KamLAND data corresponding to statistics 
of 2.6 kTy, is also shown. Since a new reactor power plant Shika-2 is planned to start operating 
in March of 2006, we have also included the contribution of the flux from this new reactor to 
the data collected in KamLAND after March 2006. The impact of the new Shika-2 reactor on 
the KamLAND sensitivity to Am^i and sin 2 #12 was studied in detail in [29] . As was shown in 
[29], with Shika-2 reactor running, the precision with which sin 2 612 can be determined in the 
KamLAND experiment does not improve considerably, while the high-LMA - low-LMA solution 
ambiguity increases. 

In the middle (right-hand) panel of Fig. 2 we show the expected allowed region when we 
combine a higher precision global solar neutrino data, including projected SNO III results, with 
prospective KamLAND data, simulated at the current low-LMA (high-LMA) best-fit point. We 
expect KamLAND to have collected approximately 1.3 kTy of data by the end of 2006, when the 

6 The total projected error in the SNO NC event rate measurement using Helium counters in the phase III of 
the experiment is expected, according to SNO, to be about 6% [39]. For the CC even rate we assume that the 
statistical error during the phase III would be approximately the same as in each of the earlier two phases, while 
the systematic error is taken to be of 4.5%, i.e., slightly smaller than the 5% reported in phases I and II of SNO. 
Thus, we assume that the total error in CC event rate measurement from all the three phases combined will be 
about 5%. 
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SNO experiment is foreseen to conclude. In Fig. 2 we present results obtained by combining the 
solar neutrino data, including SNO data with 5% and 6% total uncertainties on the measured CC 
and NC event rates, respectively, with simulated 1.3 kTy KamLAND data. We again note that 
the allowed range of Am^ diminishes remarkably with the incorporation of the KamLAND data 
simulated at the low-LMA best fit point. With the KamLAND data simulated at the high-LMA 
best fit point, the low-high LMA solution ambiguity "appears" at 95% C.L. The range of allowed 
values of sin 2 # 12 , however, in both cases remains practically unchanged with the inclusion of the 
KamLAND data. 

3 Measurement of the Solar Neutrino Oscillation Param- 
eters in the SK-Gd Reactor Experiment 

Super-Kamiokande is the world's largest running water Cerenkov detector, situated in the Kamioka 
mine in Japan. Next to SK in Kamioka is the KamLAND detector. KamLAND uses 1 kton of 
liquid scintillator to detect the V e coming from the nuclear reactors. SK receives the same reactor 
antineutrino flux as KamLAND. Since SK has ultra pure water as its detector material, it cannot 
efficiently tag the V e capture on protons. In particular, the neutron released in u e + p — ■> e + + n 
has to thermalise and then get captured to release 7, which can be then detected and the process 
tagged through delayed coincidence. The capture cross-section of thermal neutron on free protons 
in extremely small and hence SK is unable to separate the reactor v e capture events from the 
background. However, this could change if gadolinium was mixed with the SK water, as proposed 
recently in ref. [33]. The capture cross-section of thermal neutron on gadolinium is known to be 
remarkably large. In addition, neutron capture on Gd produces 8 MeV energy in photons making 
it easier for SK-Gd to detect them. Thus, SK loaded with Gadolinium could be used as a very 
big reactor anti-neutrino detector [33]. With its 22.5 kton of ultra pure water, the SK detector 
has about 1.5 x 10 33 free protons as target for the antineutrinos coming from various reactors in 
Japan. The KamLAND detector has only 3.46 x 10 31 free target protons [11]. Therefore for the 
same period of measurements, SK-Gd detector is expected to have about 43 times the statistics 
of the KamLAND experiment. The number of positron events in the SK-Gd detector is given by, 

N SK-Gd = N P J dE ms J dE u a(E u )R(E vis , E v ) £ ^ife - *e) (2) 

where E vis is the measured visible energy of the emitted positron, when the true visible energy, 
Ej is = E u — 0.80 MeV, E u being the energy of the incoming P e , a{E v ) is the u e + p — > e + + n 
reaction cross-section, Si(E u ) denotes the v e flux from the ith reactor, Lj is the distance between 
the ith reactor and Kamioka, R(E vis , E u ) is the energy resolution function of the detector, N P 
are the number of protons in the target, and Pi(P e ~> Ve) is the survival probability of the v e 
coming from the reactor %. Since the emitted positron from the antineutrino capture on proton 
will behave similarly to the elastically scattered electron by solar neutrinos, we use the energy 
resolution function for the solar neutrino detection provided by the SK collaboration [7]. The SK 
detector is located very close to KamLAND and we take the effective flux Si(E u ) and distance L; L 
for SK-Gd detector to be the same as those for KamLAND. 
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Figure 3: The 90%, 95%, 99%, 99.73% C.L. allowed regions in the Am^ — sin 2 #12 plane from an 
analysis of prospective data, obtained in 3 years of running of the SK-Gd detector with all the 
main 17 reactors running (see text for details). 



We simulate the data expected after 3 years of running of the proposed SK-Gd detector at a 
point in the low-LMA region with Am^ = 7.2 x 1CT 5 eV 2 and sin 2 #12 = 0.3; for the high-LMA 
solution the data is simulated in the point Am^ = 1.5 x 10~ 4 eV 2 and sin 2 6*12 = 0.3 1 . We 
construct a 18 bin data with a visible energy threshold of 3 MeV 8 and with bin width of 0.5 MeV. 
We wish to statistically analyze this prospective data. We define a x 2 function given by 

x 2 = J2{N^ ta - Nl heory ){al)-\N^ ata - Nf eary ) , (3) 

where N" (a = data, theory) is the number of events in the i th bin, afj is the covariant er- 
ror matrix containing the statistical and systematic errors and the sum is over all bins. The 
KamLAND experiment in their first published results of the measurement of reactor v e flux have 
reported a systematic error of 6.48% [11]. However, while KamLAND is a liquid scintillator de- 
tector, SK-Gd would be practically a water Cerenkov detector. Therefore the systematics of the 
SK-Gd experiment would be different from those of the KamLAND experiment. We would expect 
the SK-Gd systematics to be similar to those of the SK solar neutrino experiment, since even the 
energy range of the SK-Gd reactor v e experiment could be approximately the same as that of the 
SK solar neutrino experiment. We therefore assume a 5% systematic error in our analysis [7]. 

7 Unless otherwise stated, all future simulated data in this paper correspond to Am^ = 7.2 x 10~ 5 eV 2 (low- 
LMA) or Aml x = 1.5 x 1(T 4 eV 2 (high-LMA) and sin 2 6 12 = 0.3. 

8 In order to facilitate comparison with the KamLAND data, we present all results in terms of the visible energy. 
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Figure 4: The ratio of the numbers of the prospective e + — events at the SK-Gd and KamLAND 
experiments in the cases of oscillations and of absence of oscillations, with their expected la 
errorbars. The errorbars for the SK-Gd detector correspond to 3 years of data, while those for 
KamLAND are for a total statistics of 3 kTy. The upper panels are obtained assuming that the 
current 16 main reactors in Japan are operative. The bottom panels correspond to the case of 
the Shika-2 reactor operating along with the current 16 main reactors. The dashed (solid) lines 
correspond to the low-LMA (high-LMA) solution. 

Figure 3 shows the allowed regions in the Am^ — sin 2 # 12 plane, obtained using 3 years of 
thus simulated data in the SK-Gd experiment, when the true solution is in the low-LMA region. 
In obtaining this figure we have assumed that all reactors in Japan, relevant for the analysis, 
including the forthcoming Shika-2 reactor, would be running with their full power. We note that 
the size of the allowed areas diminishes significantly in both Am^ and sin 2 #12 compared to the 
currently allowed regions. In Table 1 we show the uncertainties in Am^ and sin 2 #12 after three 
years of running of SK-Gd alone. We also give the expected uncertainties if data from 5 years of 
running of the SK-Gd detector are used in the analysis. 

As Fig. 3 shows, the spurious high-LMA solution would be completely ruled out by the SK-Gd 
data, if low-LMA is the true solution. We have checked that SK-Gd can rule out the "wrong" 
solution with about 0.4 — 0.5 years of data. This is in sharp contrast with what we have obtained 
for KamLAND, where it is hard to resolve the low-LMA - high-LMA ambiguity even with 2.6 kTy 
of data once the Shika-2 reactor is switched on. To facilitate the understanding of this result we 
present in Fig. 4 the spectral distortions predicted to be observed in the SK-Gd (left panels) and 
in KamLAND (right panels) experiments. The upper panels show the ratio of the events (or event 
rates) in the cases of oscillations and of absence of oscillations, with the current main 16 reactors 
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running, while the lower panels are for the case where the Shika-2 reactor is operating along with 
the other 16 reactors. The dashed black lines are the expected spectra for the low-LMA solution, 
while the red solid lines are the spectra for the high-LMA solution. The errorbars correspond to 
the 1(7 statistical errors for 3 years of running of the SK-Gd detector and for 3 kTy (~ 7 years) 
of running of the KamLAND detector. The bin sizes used are 0.5 MeV for SK-Gd and 0.425 
MeV for KamLAND. The threshold for KamLAND is taken as Ef is = 2.6 MeV [11], while for 
SK-Gd we use a threshold of E^ s = 3.0 MeV. For the SK-Gd detector the observed spectrum 
extends to nearly 12 MeV due to the worse energy resolution of SK. Clearly, in the 16 reactor 
case the SK-Gd experiment can distinguish between the low-LMA and high-LMA spectra. The 
KamLAND experiment could also distinguish between the two solutions with a 3 kTy statistics 
(see also [26, 29, 40]) 9 . The differences between the e + — spectra in the cases of the low-LMA 
and the high-LMA solutions, measured in SK-Gd and KamLAND experiments, diminish once the 
Shika-2 reactor is switched on. The low-LMA and high-LMA solution spectra measured with the 
KamLAND detector would be almost overlapping, as can be seen in the lower right-hand panel. 
This results in the reappearance of the "wrong" high-LMA solution, as stressed in [29]. Even for 
the SK-Gd detector, the difference between the spectra corresponding to the two solutions reduces 
considerably, especially in the statistically most relevant bins. However, there are still bins at 
intermediate energies at which the two solutions can be easily distinguished. This is why the SK- 
Gd experiment could still resolve the ambiguity between the low-LMA and high-LMA solutions 
with Shika-2 reactor operating, while KamLAND cannot. 

The mixing angle is most precisely determined when the baseline corresponds to a minimum 
in the survival probability (SPMIN) in the statistically most relevant part of the positron energy 
spectrum [26]. It was pointed out in [29] that for the low-LMA solution the baseline from the 
Shika-2 reactor to Kamioka is such that the SPMIN for the Shika-2 flux occurs at E vis = AMeV . 
The Kashiwazaki-Kariwa power plant, which is the most powerful reactor complex in the world, 
is at a distance of about 160 km from Kamioka. Hence, the most significant part of the flux of u e 
from this most important source falls at Kamioka in the region of a survival probability maximum 
(SPMAX). 

The effective flux at Kamioka in absence of oscillations is defined as 

**-=4fb (4) 

For the Kashiwazaki complex it corresponds to $> KL « 7.3/z W/cm 2 , while for the Shika-2 reactor 
it is substantially smaller, §kl ~ A.lfiW/cm 2 . Thus, the contributions to the signals observed 
in SK-Gd and KamLAND due to the Kashiwazaki v e flux will be considerably larger than that 
due to the Shika-2 flux, when the contributions from all the reactors are combined. Therefore 
the resultant e + — spectrum due to the cumulative u e flux at Kamioka, produced by all relevant 
reactors in Japan, still exhibits the effect of the v e survival probability maximum. This is seen 
also in Fig. 4. 

We therefore consider a fictitious scenario in which the Kashiwazaki complex is completely 
switched off, while the Shika-2 is operating. In figure 5 we show the allowed regions obtained 

9 Let us note that using combined data from the KamLAND detector, corresponding to statistics of lkTy, and 
from the solar neutrino experiments, can permit to rule out the high-LMA solution [19, 29]. 
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Figure 5: Allowed regions in the Am^ — sin 2 9 12 plane from an analysis of prospective data from 
the SK-Gd detector collected over a period of 3 years, with the Kashiwazaki-Kariwa complex 
switched off and all other reactors, including the Shika-2 one, running. 

using the data collected by the SK-Gd detector after 3 years of operation with Kashiwazaki 
reactor complex switched off, but with the Shika-2 reactor working. We show the uncertainty on 
Am^i and sin 2 #12, determined with this set up, in Table 1. We note that instead of diminishing, 
the uncertainty in sin 2 # 12 increases. Even the spurious high-LMA solution reappears in this case. 

The predicted spectral distortions for the low-LMA and high-LMA solutions expected to be 
observed in SK-Gd and KamLAND, when the Kashiwazaki-Kariwa power plant does not operate, 
while the Shika-2 one is running at its full power, are shown in the upper panels in Fig 6. For 
comparison we also show the spectrum expected with all the main 17 reactors operating at their 
full strength. A comparison of the upper and lower left-hand panels indicates that with the 
Kashiwazaki complex switched off, the distortions in the expected spectrum for the low-LMA 
solution become smaller. We note that even though the spectra corresponding to the high-LMA 
and low-LMA solutions and measured with the SK-Gd detector are less overlapping, the difference 
in the spectral shape for the two solutions diminishes when the Kashiwazaki plant is not operating. 
Thus, it could be possible with a value of Am^ in the high-LMA zone, to describe the data 
simulated at the "true" low-LMA solution point. However, in this case a somewhat larger value of 
sin 2 9 12 is required in order to explain the larger reduction of the event rate corresponding to the 
value of Am 21 of the "true" low-LMA solution. This is what we get in Fig. 5, where the spurious 
high-LMA solution gets allowed with a higher value of sin 2 6 I 12 . 

To stress the importance of having an experimental set up in which there are no cancellations 
between SPMIN and SPMAX due reactors at very different distances, we consider a hypothetical 
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Figure 6: The ratio of the number of e + — events in the cases of oscillations and of absence of 
oscillations with their la errorbars, expected at SK-Gd (3 years of data) and KamLAND (3 kTy 
of data) experiments. The Shika-2 reactor is supposed to be operating. The upper panels are 
obtained assuming that the Kashiwazaki-Kariwa power plant is completely switched off, while 
the lower panels correspond to all the 17 reactors, including Kashiwazaki-Kariwa and Shika-2, 
working. The dashed (solid) lines correspond to the low-LMA (high-LMA) solution. 

situation where only one reactor is operating. As a concrete example we consider a scenario where 
only the Shika-2 reactor is operating. The flux from this reactor would produce a SPMIN in the 
e + — spectrum detected at SK-Gd if low-LMA was the correct solution. We show in Table 1 the 
range of allowed parameter values and the corresponding spread at 99% C.L. obtained with 7 years 
of data recorded in SK-Gd, when only Shika-2 is operating. The improvement in the precision of 
sin 2 9 12 measurement is seen to be remarkable. 



4 Combined Impact of Solar and SK-Gd Data 

In Fig. 7 we present the allowed regions obtained by a combined analysis of the current global 
solar neutrino data and prospective SK-Gd data, expected after 3 years of running of the proposed 
SK-Gd experiment. The left-hand panel shows the case of low-LMA being the true solution of the 
solar neutrino problem; we simulated the SK-Gd data at Am^ = 7.2 x 10~ 5 eV 2 and sin 2 9\i = 0.3. 
The right-hand panel illustrates the scenario we would witness if against the current trend, the 
next set of data from the KamLAND experiment would favor the high-LMA solution. The SK-Gd 
data for the right-hand panel is simulated at Am^ = 1.5 x 10 -4 eV 2 and sin 2 12 = 0.3. 
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Figure 7: The 90%, 95%, 99%, 99.73% C.L. regions in the Am 2 21 - sin 2 12 plane, allowed by the 
combined data from the solar neutrino experiments and from the SK-Gd detector after 3 years 
of running of the latter. The left-hand panel shows the allowed regions obtained in the case the 
SK-Gd data corresponded to the low-LMA solution, while the right-hand panel shows the allowed 
regions if high-LMA was the true solution. 

In the case of the low-LMA assumed to be the true solution of the solar neutrino problem, the 
inclusion of SK-Gd data from 3 years of measurements in the global solar neutrino analysis rein- 
forces the low-LMA solution. The best-fit and the range of allowed values of Am^ are determined 
almost solely by the SK-Gd data. The best-fit and the range of allowed values of sin 2 9 12 depend 
on both the solar and the SK-Gd data, though again due to its enormous statistics the SK-Gd 
data have larger impact on the sin 2 #12 determination. 

If, however, the next KamLAND results do conform to the high-LMA solution, we would en- 
counter a conflicting trend where the description of the solar neutrino data and of the KamLAND data 
would require different values of Am^. Such a situation might warrant an explanation which in- 
cludes new phenomenon along with matter induced oscillations to be responsible for the solar 
neutrino deficit (see, e.g., [41]). The data collected by SK-Gd experiment in 3 years, when com- 
bined with the global solar neutrino data, would lead to the allowed regions in the Am^ — sin 2 9 12 
plane, shown in the right-hand panel in Fig. 7. We note that if the SK-Gd data conformed to a 
point in the high-LMA zone, the low-LMA solution which is favored by the solar neutrino data, 
would be completely ruled out. This reflects the statistical power of the SK-Gd experiment. In 
a similar analysis with prospective KamLAND data only (and no SK-Gd data) corresponding to 
the high-LMA solution, we would get allowed zones in both the low-LMA and high-LMA regions 
(see, e.g., ref. [19]). 
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Figure 8: The same as in Fig. 7, but with errors on the CC and NC event rates measured at SNO 
reduced to 5% and 6%, respectively. 



In Table 1 we present the range of allowed Am^ and sin 2 9u values, determined in a combined 
analysis of the solar neutrino and SK-Gd data. In the case of the true value of Am^ lying in the 
low-LMA region, the lower limit on Am^ is less stringent in comparison with the lower limit one 
obtains from the 3 year data of SK-Gd only This is due to the solar neutrino data favoring values 
of Am^ ~6x 10~ 5 eV 2 . The upper limit on Am^ is determined solely by the SK-Gd data. The 
maximal allowed value of sin 2 # 12 diminishes with the inclusion of the solar neutrino data in the 
analysis together with the SK-Gd data. Correspondingly, the spread reduces from 19% to 15%. 
The lower limit on sin 2 # 12 is determined by the SK-Gd data. 

If the SK-Gd data is simulated at a point in the high-LMA region, Am^ has a spread of 3% 
and sin 2 9\ 2 has a spread of 12%. We note that the spread in sin 2 6*12 is larger for values of Am 2 x 
in the high-LMA region than for Am^ lying in the low-LMA region. The accuracy in the Am^ 
determination is somewhat worse in the case of the high-LMA solution compared to the case of the 
low-LMA solution because there is less spectral distortion for the Am^ values in the high-LMA 
zone, as can be seen in Fig. 4. For a high precision determination of sin 2 O12, conditions as close 
to SPMIN as possible are required [26, 29]. We find that the precision in sin 2 #12 is worse for the 
high-LMA solution since the spectrum has a smaller "dip" even for the higher energy bins, as 
compared to the "dip" in the spectrum in the case of the low-LMA solution (this is valid even if 
the Shika-2 reactor is operating, see the lower left-hand panel in Fig. 4). 

In Fig. 8 we present the the allowed regions obtained by a combined analysis of the global 
solar neutrino data with reduced errors for the SNO CC (5%) and SNO NC (6%) and prospective 
SK-Gd data, expected after 3 years of running of the proposed SK-Gd experiment. In Table 1 we 
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present the corresponding allowed range for Am^ and sin 2 612 and the spread at 99% C.L. 

5 Conclusions 

In this paper we have studied the prospects of high precision determination of the solar neutrino 
oscillation parameters, Am^ and sin 2 9 12 , by investigating the reactor u e oscillations with the 
recently proposed gadolinium loaded Super-Kamiokande (SK-Gd) experiment [33]. Doping the 
SK detector with Gadolinium by dissolving 0.2% gadolinium trichloride in the water [33] would 
make it possible to detect the neutrons, released in the v e capture on protons, with a relatively 
high efficiency This would allow the SK experiment, in particular, to detect the reactor v e coming 
from the numerous powerful nuclear reactors located in Japan, thus transforming SK into the 
largest long baseline reactor v e detector with an event rate approximately 43 times higher than 
that observed in the KamLAND experiment. 

Working in the framework of two-neutrino mixing, we have presented first an overview of the 
currently allowed ranges of values of the solar neutrino oscillation parameters Am^ and sin 2 9\2 
by the existing global solar neutrino and KamLAND data (Fig. 1 and Table 1), as well as the 
regions in the parameter space which are expected to be allowed after taking into account future 
higher statistics and lower systematic data from SNO and from KamLAND experiments (Fig. 
2). The SK-Gd detector is considered in Section 3, where we discuss the assumptions made to 
simulate the data in this proposed modification of the SK detector, as well as the procedure used 
to statistically analyze the simulated data. 

The results of our analysis show that the SK-Gd experiment has a remarkable potential in re- 
ducing the uncertainties in the values of Am^ and sin 2 9 12 (Figs. 3 and 7 and Table 1). Combining 
the SK-Gd data, taken over a period of 3 years, with the data from the solar neutrino experiments, 
will allow to determine Am^ with a 3% error at 99% C.L., the error being by a factor of 2 smaller 
if one uses only the SK-Gd data. With only ~ 0.5 years of data, the SK-Gd experiment would 
rule out completely the spurious high-LMA solution if low-LMA is the true solution. In contrast, 
the high-LMA solution could be ruled out if one uses the data from the solar neutrino experi- 
ments and data from the KamLAND detector, corresponding to lkTy of statistics. If the future 
KamLAND data would conform to a point in the high-LMA region, both the high-LMA and the 
low-LMA solutions would be allowed at 90% C.L. [19]. In this case the high-LMA - low-LMA 
solution ambiguity can be resolved completely by the SK-Gd experiment. In what regards the 
solar neutrino mixing angle, using the data from the solar neutrino experiments and the SK-Gd 
3 year data will allow to determine sin 2 #12 with a ~ 15% error at 99% C.L. One would have the 
same uncertainty in sin 2 12 using data only from the SK-Gd experiment, corresponding to 5 years 
of measurements. Thus, the SK-Gd experiment could provide one of the most precise (if not the 
most precise) determination of the solar neutrino oscillation parameters Am 2 x and sin 2 6*12. 

Addendum after the 766.3 kTy results from KamLAND 

The KamLAND collaboration has published data corresponding to a statistics of 766.3 kTy [42]. 
The new results from KamLAND have corroborated their first results from 2002 and have for 
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Figure 9: Same as Fig. 3 but with the SK-Gd data simulated at the new global best-fit point, 
Am^ = 8.3 x 10~ 5 eV 2 and sin 2 12 = 0.27, marked on the figure by a star. 

the first time provided unambigous evidence for distortion in the measured reactor antineutrino 
spectrum. These new results have reinforced the low-LMA solution of the solar neutrino problem 
and the combined solar and KamLAND data now rule out the high-LMA solution at the 4<r level 
[43]. The best-fit values of the solar neutrino oscillation parameters obtained from the combined 
analysis of solar and KamLAND data are Am 21 = 8.3 x eV 2 and sin 2 9yi = 0.27 and their 
corresponding 99% CL. allowed range of values are given as 

Am 2 x = (7.3 - 9.4) x 10~ 5 eV 2 ; spread = 13% (5) 
sin 2 ^i 2 = (0.22 - 0.36); spread = 24% (6) 

We present in this section the allowed regions in the Am 21 — sin 2 9 12 parameter space expected 
from the Sk-Gd experiment when the "data" is simulated at the new global best-fit point Am 21 = 
8.3 x 10 -5 eV 2 and sin 2 # 12 = 0.27. The Fig. 9 shows the allowed regions for 3 years statistics in 
SK-Gd. The corresponding allowed range of Am 21 and sin 2 9 U and their spread, expected at 99% 
C.L. is given as 

Am 2 x = (8.01 - 8.61) x 10" 5 eV 2 ; spread = 3.6% (7) 
sin 2 n = (0.22 - 0.34); spread = 21% (8) 

Fig. 10 shows the allowed regions for 5 years statistics in SK-Gd with the spectrum simulated 
at the new global best-fit. The 99% C.L. allowed range and spread in Am 21 and sin 2 9 12 are given 
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Figure 10: Same as Fig. 9 but for 5 years statistics. 

as 

Am 2 21 = (8.07 - 8.53) x 10~ 5 eV 2 ; spread = 2.8% (9) 
sin 2 # 12 = (0.22 - 0.32); spread = 18% (10) 

We note that the precision on Am^ and sin 2 # 12 measurement in SK-Gd for a given statistics 
remains roughly the same for the old and new global best-fit points. 
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